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Abstract. The minimum and maximum mass of protoneutron stars and neutron stars are investigated. The hot 
dense matter is described by relativistic (including hyperons) and non-relativistic equations of state. We show 
that the minimum mass (~0.88- 1.28 Mq) of a neutron star is determined by the earliest stage of its evolution 
and is nearly unaffected by the presence of hyperons. The maximum mass of a neutron star is limited by the 
protoneutron star or hot neutron star stage. Further we find that the delayed collapse of a neutron star into a 
black hole during deleptonization is not only possible for equations of state with softening components, as for 
instance, hyperons, meson condensates etc., but also for neutron stars with a pure nucleonic-leptonic equation of 
state. 
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1. Introduction 

It is expected that a neutron star (NS) is born as a result 
of the gravitational collapse of the iron cor e of a massive 
(M > 8M ) evolved star (e.g. |Bethe 1990| ). Shortly after 
core bounce (some 10 ms) a hot, lepton rich NS, called 
protoneutron star (PNS), is formed. This PNS evolves in 
some minutes into a cold, lepton pure NS due to the loss 
of neutrinos (e.g. Burrows & Lattimer 1986; Keil & Janka 



|1995| ; |Pons et al. 1999| ). During these first minutes the PNS 
can collapse delayed to a black hole (BH) if its initial mass 
is high enough, either due to the loss of neutrinos or due 
to post bounce accretion. 

The aim of this work is to study the minimum and the 
maximum mass of NSs and the possibility of BH forma- 



thcsc investigations start with models of PNSs about 1 - 
3 s after core bounce (see Sect. || for more details). In this 
paper the earliest stage of a PNS is taken additionally into 
account, which is reached about 100 ms after core bounce, 
to calculate the maximum mass of a NS (see Sect. | for 
more details). 

2. Inside a protoneutron star 

As pointed out in the Introduction a special topic of this 
investigation is the more detailed incorporation of the ear- 
liest stage of the PNS, about 0.1-1 s after core bounce 



(Burrows & Lattimer 1986|; 


Burrows et al. 1995; 


Keil et al. 


199G: 


Pons et al. 1999 


). This early type PNS is charac- 



tion during the deleptonization period. Recently Gondek 
et al. (1997|, |1998| ), |Goussard et al. (1998| ) and |Strobel 



et al. (1999a) have calculated limits on the minimum mass 
of NSs, using equations of state (EOSs) including nucle- 
ons (n, p) and leptons. In this paper refinements to these 
approaches are given by calculating PNSs for a large sam- 
ple of EOSs using a non-relativistic Thomas-Fer mi model 
([Myers fc Swiatecki 1996] ; jStrobcl et al. 1999b|Ja|) and a 



tcrized by a hot shocked envelope with an entropy per 
baryon of s ~4-6 (in units of the Boltzmann constant, 
fee) for baryon number densities n < n onv , an unshocked 
core with s ~l-2 for densities n > n corc , and a transi- 



tion region between these densities (Burrows et al. 1995; 
Keil et al. 1996| ). For the description of this early stage 



of the PNS we use, in accordance with the investigations 
of [Burrows & Lattimer (1988), Burrows et al. (1995| ) and 



relat ivistic Hartree model ( perot fc Walecka 1986[ p chafer 



1997 ) including hyperons. Limits on the maxi mum mass 



1995; 


Bombaci 1996 


Gondek et al. 199? 



of NSs were studied by numerous aut hors (e.g. Takatsuka 
Ellis et al. 1996| ; |Prakash et al. 1997| ; 



using different kinds of EOSs. All 
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Strobel, 



Keil et al. (1996), baryon densities in the range of 0.002 
0.02 fm~ 3 for n cnv and 0.1 fm -3 for n Gore . The entropy per 
baryon is chosen between 5 - 6 in the envelope and be- 
tween 1 - 2 in the core, respectively. The investigation is 
performed for the early type PNS models with trapped 
neutrinos using a constant lepton number (Yj = 0.4) for 
densities above n = 6 x 10~ 4 fm~ 3 . 

We assume that post bounce accretion onto the pro- 
toneutron star stops within the first 500 ms after core 
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Table 1. Properties of cold symmetric nuclear matter. 
The entries are: energy per baryon, u; saturation density, 
no; incompressibility, K^; symmetry energy, J; effective 
nucleon mass, m* jm. 



u no Koo J m*/m 

[MeV] [fm- 3 ] [MeV] [MeV] 



Table 2. Coupling constants for the relativistic Hartree 
models. The coupling constants are given for the following 
meson masses: m a — 550 MeV, m u — 783 MeV and 
nip — 770 MeV (the coupling constants are converted 
according to these meson masses). 



ff?/47r gt/A-K g 2 p /4ir 10 3 b 10 3 c 



TF 


-16.24 


0.161 


GM1 


-16.30 


0.153 


GM3 


-16.30 


0.153 


NL1 


-16.42 


0.152 



234 


32.7 


0.87 


300 


32.5 


0.70 


240 


32.5 


0.78 


212 


43.5 


0.57 



GM1 7.288 8.959 
GM3 6.139 6.041 
NL1 10.2099 13.6108 



5.346 2.947 -1.070 
5.807 8.659 -2.421 
8.0260 2.4578 -3.4334 



bounce. The ammount of accreted matter, and the time 
at which the accretion through the shock front stops are 
still open questions. The amount of accreted matter after 
core bounce ranges from 0.5 Mq during the first 10 ms or 
so to 0.001 to 0.15 Mq during the following 100 ms or so, 
but in most calculations significant accretion stops about 



topic see e.g. 


Burrows & Lattimcr 1988 


Chevalier 1989; 


Herant et al. 1992; Burrows et al. 1995; 


Janka & Miiller 


1996; 


Mezzacappa et al. 1998 




Fryer & Heger 2000 


)• 



About 1 - 3 s after core bounce the PNS has a nearly 
constant entropy per baryon (s ~ 2) and an approxi- 
mately constant lepton number (Y\ ~ 0.4) for densities 
n > 6 x 10~ 4 fm~ 3 . The next stage in the lifetime of 
a NS is reached after about 10-30 s, where the hot NS 
(s ~ 2) is now transparent to neutrinos. During the fol- 
lowing minutes this hot NS evolves into a cold NS. For a 
more detailed description of the evolution of PNSs see, for 
instance, Prakash et al. (1997) and Strobel et al. (1999a| ). 

The EOSs used in this paper for the description of 
PNSs are: (i) A non-relativistic Thomas-Fermi model 



(TF) for finite temperatures developed by Strobel et al 
(1999a[|]). (ii) A relativistic Hartree model for finite tem- 
peratures developed by Schafcr (1997j ). The Hartree model 
includes the following hyperons and nucleonic isobars: 
A, E~, E°, £ + , S~, S°, A - , A , A + , A ++ (the Delta reso- 
nances are not included in the NL1 model). 

The main properties of cold symmetric nuclear matter 
of the differe nt models are listed in Table [T| (s ec Strobel 



et al. 1999a; Glcndcnning & Moszkowski 1991; Reinhard 
1989[)~ The coupling constants for the Hartree models 
(GM1, GM3, NL1) are listed in Table |. We take the 
nucleon-hyperon coupling strength equal to the nucleon- 
nucleon coupling strength for simplicity. This restriction 
will not change the results for the minimum mass of a 
NS, as we will show later. The influence on the maxi- 
mum mass will be strong, but for this point we refer to 
Huber et al. (1998 ), where this question is studied for a 
large spectrum of different cmtpliag constants and Hiffpr- 

ent nnitlnnn-liypprnn r-nnpling strpncrttis fnr i-nlH NSs Ttin 

EOSs for the different models with nucleonic-leptonic mat- 
ter (TF, GM1, GM3 and NL1) and nuclear matter includ- 




Fig. 1. Pressure versus baryon number density in the high 
density region for different EOSs. The solid curves corre- 
spond (from the bottom to the top) to the cold EOS, the 
s = l,Yi = 0.4 EOS and the s — 2,Y\ = 0.4 EOS of the 
TF model. The long dashed curves correspond to the cold 
EOS (lower) and the s = 1, Y } = 0.4 EOS (upper) of the 
NL1 model without hyperons (dot-dashed lines, hyperons 
included) . 



ing hyperons (GMlnyp, GM3H yp and NLlnyp) are shown 
in Fig. [l] and Fig. |[ The EOSs for baryon number densi- 
ties below n = 0.1 fm -3 are taken from the investigation 



of Strobel et al. (1999a 



The NL1 coupling constants lead to the stiffest EOSs 
at high densities, since the high density behaviour is dom- 
inated by the strength of the w-coupling. However the in- 
compressibility is the smallest. In this respect one has to 
keep in mind that the incompressibility determines the 
stiffness only near saturation and for higher densities the 
second and third derivatives of the incompressibility are 



important for the stiffness (see for instance Ramschiitz 



et al. 199C) 
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Table 3. EOSs used in this paper and resulting minimum and maximum masses. The entries are: time after core 
bounce, i; entropy per baryon in the envelope, s env ; entropy per baryon in the core, Score! maximum baryon number 
density of the envelope correlated with the entropy per baryon in the envelope, n en v(senv); minimum baryon number 
density of the core correlated with the entropy per baryon in the core, n cor e(score); minimum baryonic mass of the 
PNS, M™ n ; resulting minimum gravitational mass of the cold NS, M™ m (T = 0); maximum baryonic mass of the 
PNS or NS, M^ ax ; maximum gravitational mass of the cold NS determined by the PNS or hot NS, Mg lax (T = 0); 
possibility of the formation of a black hole during deleptonization, BH. 



Model 


t 






^env (Senv) 


'"'core (^corc ) 


MS 1 '" 


MJ? in (T - 0) 


ivij3 




BH 




L°J 






ffm~ 3 l 
[nil j 


ffm" 3 l 
[mi j 


\Mr-A 

i ©J 




V 1V1 Q\ 


I 0J 




TF 


0.1 - 1 


6.0 


2.0 


0.002 


0.1 


1.338 


1.284 


2.366 




yes 




0.1 - 1 


6.0 


1.0 


0.002 


0.1 


0.949 


0.878 


2.371 




yes 




0.1 - 1 


6.0 


1.0 


0.005 


0.1 


1.117 


1.021 


2.372 




vos 
j *■ ° 




0.1 - 1 


5.0 


1.0 


0.01 


0.1 


1.020 


0.940 


2.372 




y eb 




0.1 - 1 


■3.0 


1.0 


0.02 


0.1 


1.229 


1.114 


2.374 




VPS 

y co 




1 - 3 


2.0 


2.0 










2.364 


1.969 






10 - 30 


2.0 


2.0 










2.449 








oc 


0.0 


0.0 










2.417 






GM1 


0.1 - 1 


6.0 


1.0 


0.005 


0.1 


1.255 


1.206 


2.629 




ves 




0.1 - 1 


5.0 


1.0 


0.005 


0.1 


1.019 


0.988 


2.628 




ves 




0.1 - 1 


■3.0 


1.0 


0.01 


0.1 


1.140 


1.100 


2.629 




VPS 

y co 




1 - 3 


2.0 


2.0 










2.613 


2.251 






10 - 30 


2.0 


2.0 










2.670 








OC 


0.0 


0.0 










2.706 






GMlHvn 


0.1 - 1 


6.0 


1.0 


0.005 


0.1 


1.255 


1.206 


2.462 




ves 




0.1 - 1 


5.0 


1.0 


0.005 


0.1 


1.019 


0.987 


2.461 




ves 




0.1 - 1 


■3.0 


1.0 


0.01 


0.1 


1.140 


1.100 


2.462 




VPS 




1-3 


2.0 


2.0 










2.401 




ves 




10 - 30 


2.0 


2.0 










2.304 


2.015 






oo 


0.0 


0.0 










2.353 






GM3 


0.1 - 1 


6.0 


1.0 


0.005 


0.1 


1.255 


1.205 


2.204 


1.939 


no 




0.1 - 1 


•3.0 


1.0 


0.005 


0.1 


1.039 


1.005 


2.203 


1.938 


no 




0.1 - 1 


5.0 


1.0 


0.01 


0.1 


1.145 


1.104 


2.204 


1.939 


no 




1 - 3 


2.0 


2.0 










2.212 








10 - 30 


2.0 


2.0 










2.249 








oo 


0.0 


0.0 










2.247 






GM3 H yp 


0.1 - 1 


6.0 


1.0 


0.005 


0.1 


1.255 


1.205 


2.047 




yes 




0.1 - 1 


5.0 


1.0 


0.005 


0.1 


1.039 


1.005 


2.046 




yes 




0.1 - 1 


5.0 


1.0 


0.01 


0.1 


1.145 


1.104 


2.047 




yes 




1 - 3 


2.0 


2.0 










2.015 




yes 




10 - 30 


2.0 


2.0 










1.898 


1.699 






oo 


0.0 


0.0 










1.911 






NL1 


0.1 - 1 


5.0 


1.0 


0.005 


0.1 


1.175 


1.111 


3.218 




yes 




1 - 3 


2.0 


2.0 










3.162 


2.663 


(yes) 




10 - 30 


2.0 


2.0 










3.202 


(2.689) 






oo 


0.0 


0.0 










3.297 






NLlHyp 


0.1 - 1 


5.0 


1.0 


0.005 


0.1 


1.169 


1.104 


3.017 




yes 




1 - 3 


2.0 


2.0 










2.942 




yes 




10 - 30 


2.0 


2.0 










2.871 


2.399 






oc 


0.0 


0.0 










2.962 







3. The minimum mass of neutron stars 

The minimum and maximum mass of non-rotating NSs 
for the different EOSs and different choices of the tran- 
sition region between the hot shocked envelope and the 
unshocked core of the PNS are given in Table [| It turned 
out that the minimum mass of a NS is determined by 
the earliest stage of the lifetime of a PNS. The values for 



the minimum gravitational (baryonic) mass, Mq (Mb), 
lie in the range of 0.878- 1.284 Mq (0.949- 1.338 Mq) for 
the different models (see Table y). They are larger by a 
factor of ten than t he minimum mass of a cold NS (e.g. 
trobel et al. 1997) emerging from a calculation for the 



I 



minimum mass of the cold NS based only on the solu- 
tion of the Tolman-Oppenheimer-Volkoff equation with 
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Fig. 2. Pressure versus baryon number density in the high 
density region for different EOSs. The dotted curves cor- 
respond to the cold EOS (lower) and the s = 1,Y\ = 0.4 
EOS (upper) of the GM1 model without hyperons (long 
dashed lines, hyperons included). The dot-dashed curves 
correspond to the cold EOS (lower) and the s = 1,1] = 
0.4 EOS (upper) of the GM3 model without hyperons 
(short dashed lines, hyperons included). 



the corresponding EOS for cold NS matter. This shift due 
to the p roperties of early typ e P NSs was also recentl y 
found by|Goussard et al. (1998b and |Strobel et al. (1999a|) . 
Gondek et al. (1997 , 1998] ) found an even smaller value 
for the minimum mass (Mq ~ 0.6 Mq). The reason for 
this is that they use models for PNSs about 1-3 s af- 



ter core bounce, where t he minimum mass is lower, see 
also Strobel et al. (1999a ). Our models with hyperons in- 



cluded do not change the values for the minimum mass 
significantly (AM G < O.OO7M ). This result could be 
expected since the maximum baryon number density in- 
side the core of such a NS is smaller than two times the 
nuclear matter density and hence the hyperon fraction is 
rather small. 

Let us now turn to the question, why there is a mini- 
mum mass in this early stage of the evolution. A star be- 
comes dynamically unstable if its mean adiabatic index, 
f , lies below 4/3. In the case of a cold NS this is caused 
by the neutron drip at a density of n ~ 2.6 x 10~ 4 fm~ 3 . 
For the earliest stage of the PNS the drop below T = 4/3 
is caused by the entropy drop between the hot shocked 
envelope and the unshocked core of the PNS. For a com- 
plete description of the stability criteria of neutron stars 



see Shapiro fc Teukolsky (1983) 



4. The maximum mass of neutron stars and the 
possibility of black hole formation during 
deleptonization 

Quite interesting is the result for the maximum gravita- 
tional mass for NSs, for which we obtain values in the 



2.8 



2.75 



2.7 



2.65 



2.6 




3.1 



3.15 



3.2 



3.25 



t = 0.1-ls 
t = 1-3 s 
t = 10-30 s 
t = minutes 



3.3 



3.35 



Fig. 3. Maximum gravitatinal versus maximum baryonic 
mass for different evolutionary stages of the NL1 EOSs. 
The two lines denote the mass window for the delayed 
collapse. 



range between 1.699- 2.663 Mq (see Table |) for the dif- 
ferent EOSs used in this paper, with the possibility of 
the formation of a BH during deleptonization for EOSs 
with a pure nucleonic-leptonic composition and EOSs 
which include hyperons. For instance, for the NL1 EOS, 
the maximum baryonic mass for the first milliseconds, 
Afg ax = 3.218 Mq, is larger than the allowed maximum 
baryonic mass of 3.162 Mq after 1 - 3 s, which implies the 
possibility of a BH formation (see Table || and Fig. ||). 
This value also sets the boundary of 2.663 Mq for the 
maximum gravitational mass for the cold NS. The maxi- 
mum baryonic mass of the first stage is even larger than 
maximum value 10-30 s after core bounce fore the NL1 
EOS. This means that the delayed collapse is also possi- 
ble for the totally deleptonized hot NS of the NL1 EOS. It 
should be mentioned in this context that the exact entropy 
profile of the hot shocked envelope of an early type PNS 
plays a minor role in determining the maximum baryonic 
mass, since a PNS with a constant entropy per baryon of 
s = 1 throughout the whole PNS and a lepton number 
of Y\ = 0.4 leads to nearly the same maximum bary- 
onic mass (e.g. for the NL1 model: Mf ax = 3.218 Mq, 
0.1-1.0 s after core bounce; Mf ax = 3.219 Mq, for an 
entropy per baryon of s = 1 throughout the whole PNS 
and Y\ = 0.4). This means that the maximum mass de- 
pents nearly complete on the properties of the core of the 
PNS0. 

Similar behaviour is found for the other EOSs with 
exception of the nucleonic-leptonic EOS GM3, where the 
allowed maximum mass after 1-3 s exeeds the maxi- 
mum masses of the first period. The reason why this 
possibility of a BH formation was dismissed so far for a 



1 Note: Other properties of an early type PNS (e.g. the ra- 
dius) are strongly affected by the entropy profile. 
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Fig. 4. Baryonic mass difference between the early type 
PNS stage (0.1-1 s after core bounce) and the late type 
PNS stage (1-3 s after core bounce) versus resulting 
maximum baryonic mass for the different pure nucleonic- 
leptonic EOSs. 



pure nuclconic-lcptonic composition (e.g. Takatsuka 1995; 
Bombaci 199(| [Ellis ct al. 1996| ; [Prakash et al. 1997| ; 
Gondck ct al. 1998Q was that these authors start their 
calculations about 1-3 s after core bounce, since the 
following stages allow only larger baryonic masses. One 
should mention in this context that for the nucleonic- 
leptonic EOSs the differences in the maximum masses 
in the first seconds are more than three times smaller^] 
(< 0.056 A/q) than for the softer EOSs with hyperons 
(0.146-0.158 Mq). Due to this, for the EOSs including 
hyperons the possibility of a BH formation shows up even 
by starting the calculation 1 - 3 s after core bounce (see 
Table ||). It turned out, that the mass window for the de- 
layed collapse of the pure nucleonic-leptonic EOSs is an 
inceasing function of the maximum possible mass of the 
different models (see Fig. |Z]). The EOSs including hyper- 
ons did not show a similar behaviour, the mass window is 
nearly constant for all used EOSs. The minimum value of 
the maximum baryonic mass for the models including hy- 
perons is reached by the hot NS without trapped neutrinos 
(10-30 s after core bounce). This, as already mentioned, 
can lead to a delayed collapse of the NS during delep- 
toniz ation if the initial mass of the PNS is high enough 
(e.g. Baumgarte et al. 1996; Pons et al. 1999| ). 

It is an open question how large the maximum mass 
of a NS really is. It is clear that an EOS should allow NS 
masses larger than 1.4 Mq since measurements of pulsars 
in binary systems show up values around this mass (e.g. 
Thorsett & Chakrabarty 1999). If quasi periodic oscilla- 
tions (e.g. van der Klis 2000) are taken into account, the 
possible maximum mass of an EOS has to be larger than 
1.8 - 2.0 M© (e.g. [Schaab fc Weigel 1999| ). This whould 



lead to the conclusion that only the GM1 and the NL1 
EOSs (with and without hyperons) are stiff enough to al- 
low NSs with a larger mass. 

5. Discussion and conclusion 

In this paper we presented a calculation of the minimum 
and maximum mass of a NS taking the early evolution, 
shortly after core bounce, into account. We found that 
the minimum mass of a NS is limited by the earliest stage 
of its evolution (~ 0.1-1 s after core bounce). Therefore 
the minimum gravitational mass has values between 0.878 
and 1.284 Mq. These values li e in the same range as 
thos e found in recent wo rks by Goussard et al. (1998 ) 
and Strobel et al. (1999a). Additionally we have exam- 



ined the influence of hyperons on the minimum mass. As 
expected it turned out that the minimum mass is nearly 
unaffected by the inclusion of hyperons. Another possibil- 
ity of NS pro duction is the accretion ind uced collapse of a 
white dwarf ( Canal fc Schatzman 1976|), but the numeri- 
cal simulation of Woosley & Baron (1992) showed that the 
resulting NS will have a mass of ~ 1.4 Mq since the white 
dwarf collapses without significant mass loss after reach- 
ing the maximum mass of a white dwarf. From the results 
reported it is obvious that NSs, which are born in Type-II 
supernovae or due to accretion induced collapse of a white 
dwarf, have a minimum mass about ten times larger than 
the possible minimum mass of a cold NS (~ O.IMq). 

With respect to the maximum mass of a NS calculated 
by solving the Tolman-Oppenheimer-Volkoff equation for 
NS matter the outcome is rather uncertain, since the EOS 
in the high density region rests strongly on theoretical ex- 
trapolations. The underlying model, for instance, different 
many-body approximations and/or inclusion of more mas- 
sive baryons, quarks, condensates etc., can rende r the re- 
sults for the maximum masses significantly (e.g. Prakash 
|ct al. 1997 ; Huber et al. 1998 ). However, as a result of 



our work it turned out, that the maximum mass of a NS 
is limited by the first ~ 30 seconds of its evolution. The 
maximum baryonic mass of a NS, constructed with pure 
nucleonic-leptonic EOSs and EOSs including hyperons, is 
always reached during this early epoch of the evolution. 

Finally we found that BH formation during delep- 
tonization is also possible for NSs with pure nucleonic- 
leptonic EOSs if the initial mass is high enough. This prop- 
erty of NSs was expected till now to be possible only for 
EOSs including softening ingredients, like hyperons, me- 
son condensates or a quark-hadron phase transition (e.g. 
Brown fc Bet he 1994| ; |Glendenning 1995|). Ther e are nu- 



merous papers on this topic (e.g. Takatsuka 1995; Bombac 



1996; Ellis et al. 1996; Prakash et al. 1997; Gondek et al 



1998), but all these investigations started their calcula- 



For that reason the effect was overlooked in strobel et al 



(1999a) 



tions about 1 - 3 s after core bounce. In this paper we 
included the earliest stage and found that the maximum 
baryonic mass of the early type PNSs (~ 0. 1 - 1 s after core 
bounce) could be up to 0.056 Mq larger than that of the 
late type PNS (~l-3 s after core bounce). Due to this 
result it seems possible that a NS with a pure nucleonic- 
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leptonic EOS could also undergo delayed collapse into a 
BH during the first seconds of its evolution. This question 
could only be finally decided in time evolution calculations 
like this including softening ingredients (e.g. Baumgarte 
et al. 996), since we have used approximate values for 



the entropy per baryon and the lepton number at fixed 



times. Post bounce accretion (e.g. Chevalier 1989 
et al 



Hcrant 



992) could lead also to a delayed collapse, but it 



is uncertain till now how much matter falls through the 
shock front within the first second after core bounce, so 
that if a cut off of a neutrino signal after a few seconds 
is detected in the future, the two ways of BH formation 
could possibly be indistinguishable. 
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